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Abstract 

An interesting issue in fluid dynamics is represented by the possible existence of inverse kinetic 
theories (IKT) which are able to deliver, in a suitable sense, the complete set of fluid equations 
which are associated to a prescribed fluid. From the mathematical viewpoint this involves the 
formal description of a fluid by means of a classical dynamical system which advances in time the 
relevant fluid fields. The possibility of defining an IKT for the 3D incompressible Navier-Stokes 
equations (INSE), recently investigated (Ellero et al, 2004-2007) raises the interesting question 
whether the theory can be applied also to thermofluids, in such a way to satisfy also the second 
principle of thermodynamics. The goal of this paper is to prove that such a generalization is 
actually possible, by means of a suitable extended phase-space formulation. We consider, as a 
reference test, the case of non-isentropic incompressible thermofluids, whose dynamics is described 
by the Fourier and the incompressible Navier-Stokes equations, the latter subject to the conditions 
of validity of the Boussinesq approximation. 

PACS numbers: 52.25.Dg,47.10.ad,05.70.Ln 
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I. INTRODUCTION 



A remarkable aspect of fluid dynamics is related to the construction of inverse kinetic 
theories (IKT) for hydrodynamic equations in which the fluid fields are identified with suit- 
able moments of an appropriate kinetic probability distribution. Recently the topic has 



been the subject o 
equations (INSE) ^ |3 



theoretical investigations on the incompressible Navier-Stokes (N-S) 



The importance of the IKT-approach goes beyond the aca- 
demic interest. In fact, fluid equations represent usually a mixture of hyperbolic and elliptic 
pde's, which are extremely hard to study both analytically and numerically. As such, their 
investigation represents a challenge both for mathematical analysis and for computational 
fluid dynamics. For this reason in the past alternative approaches, based on asymptotic 
kinetic theories, have been devised which permit to advance in time the fluid fields, to be 
determined in terms of suitable moments of an appropriate kinetic distribution function. An 
example is provided by kinetic theories for incompressible fluids which adopt the so-called 



Lattice-Boltzmann approach [8|] (see also related discussion in Refs. jlO|, lUl]). These meth- 
ods, which approximate the fluid equations only in an asymptotic sense, are based on the 
introduction of suitably modified (fluid) equations which permit to advance in time the fluid 
fields only in an approximate sense. In particular, typically, their modified fluid equations 
actually describe weakly-compressible fluids. The discovery of IKT [1] provides, however, 
a new starting point for the theoretical and numerical investigation of hydrodynamic equa- 
tions, since it does not require any modification of the exact fluid equations, in particular it 
holds for strong solutions, and permits to advance in time exactly the fluid fields by means 
of a suitable kinetic distribution function /(x,t). Here x is the state vector x = (r 1 ,v 1 ), 
where respectively ri and vi denote the corresponding "configuration" and "velocity" vec- 
tors, and T is the phase-space spanned by x. In the sequel we shall assume that T is an 
extended phase-space, i.e., it has a dimension 2n with n > 3. This is achieved introducing a 
phase-space classical dynamical system 

x x(t) = T ti i o x , (1) 
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which uniquely advances in time the fluid fields by means of an appropriate evolution oper- 
ator T t;to [J, |5|. This is assumed to be generated by a suitably smooth vector field X(x, t), 

|x = X(x, t ), (2) 
x(t G ) = x D , (3) 

Therefore, introducing the corresponding microscopic distribution function /(x,t), it fulfills 
necessarily in T the differential Liouville equation 

L/(x,t) = 0, (4) 

where L denotes the Liouville streaming operator Lf = J^/ + ■ {X(x, t)f}. This equation, 
which may be interpreted as a Vlasov-type inverse kinetic equation (IKE), can in principle 
be defined in such a way to satisfy appropriate constraint equations. In particular, thanks to 
the arbitrariness of the dynamical systems (CQ) [i.e., the arbitrariness of X(x, t)}, the velocity 
moments of /(x,i) might be identified - in principle - so that suitable velocity moments (of /) 
coincide with t 
as in Refs. 

fluid mass density, i.e., there results p = f dvi/(x,i). 

An interesting issue is whether the theory can be applied also to thermofluids. Such a 



he relevant fluid fields characterizing a prescribed classical fluid. For example, 
one can impose that the first velocity-moment coincides with the 



generalization, as shown in an accompanying paper [1 Oj j . is actually non- unique. In partic- 
ular, the goal of this paper is to prove that an IKT can be achieved by means of a suitable 
extended phase-space formulation. We consider, as a reference problem, the case of incom- 
pressible thermofluids subject, for greater generality, to the condition of non-isentropic flow. 
For definiteness, we shall assume that the relevant fluids {p = p Q > 0, V,p > 0, T > 0, St} , 
i.e., respectively the (constant) mass density, fluid velocity, pressure, temperature and en- 
tropy describing the fluid, are defined in an appropriate existence domain. In particular 
if Q is an open connected subset of R 3 (denoted as configuration domain; with prescribed 
fixed boundary Sfl and closure Q defined as the set where the mass density is a constant 
p > 0) and I a finite time interval I =]to, ti[(with closure I— [i^ti]), we assume that the 
fluid fields {p, V, T} are continuous in Q x J, satisfy suitable initial and boundary condi- 
tions respectively at t = t and on SQ, while in the open set fix I they satisfy the so-called 
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non-isentropic and incompressible Navier- Stokes- Fourier equations (INSFE), i.e., 



V- V = 0, 

J^V + V • VV + — [Vj9 - f] - z/V 2 V = 0, 




(5) 
(6) 
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1 



+ 



PoCp 



(7) 



(8) 



Here the notation is standard. Thus, Eg. ([3]) denotes the so-called isochoricity condition, 
while Eq.(j6]) is the Navier- Stokes equation in the Boussinesq approximation. Hence, in such 
a case the force density f reads f = p Q g (1 — k p T) + fi, where the first term represents the 
(temperature-dependent) gravitational force density, while the second one (fi) the action of 
a possible non-gravitational externally-produced force. Moreover: 

• Eq.(J7J) is the Fourier equation for the temperature T, to be assumed strictly positive 
in Q x /, with J the quantity of heat generated by external sources per unit volume 
and unit time (for example, Joule heating). Thus for an isolated fluid there results by 
definition J = in Q x I. 

• Eq.([8]) defines the so-called 2nd principle for the thermodynamic entropy St- For its 
validity in the sequel we shall assume that there results either everywhere in Q x J, 
J = (thermally-isolated thermofluid) or 



(externally heated thermofluid) . 

• In these equations g, k p , u, \ an d c p are all real constants which denote respectively the 
local acceleration of gravity, the density thermal-dilatation coefficient, the kinematic 
viscosity, the thermometric conductivity and the specific heat at constant pressure. 
Thus, by taking the divergence of the N-S equation (JHJ), there it follows the Poisson 
equation for the fluid pressure p which reads 




(9) 



V 2 p = -p V ■ (V • VV) + V ■ f , 



(10) 



with p to be assumed non negative and bounded in f2 x 7; 
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Finally, it is assumed that: 

• Eqs. ©-([T]) satify a suitable initial-boundary value problem (INSFE problem) and 
that a smooth (strong) solution exists for the fluid fields {p = p Q > o, V, p > 0, T > 0} ; 

• the entropy functional St can be defined so that it satisfies the 2nd principle. 



II. EXTENDED PHASE-SPACE IKT FOR INCOMPRESSIBLE THERMOFLU- 
IDS 



Here we intend to show that an IKT for INSFE can be reached by introducing of a 
suitable extended phase-space formulation, based on a generalization of the IKT developed 



previously for the incompressible Navier-Stokes equations |2 
us introduce the notations 



m. 



For definiteness, let 



n = (r,0) , 

vi = (v,w), (11) 

X = {v 1) F 1 (x,t)} ) 

Fx(x,t) = {F(x,t),tf(x,t)} 

where the vectors r and v span, respectively, the whole configuration domain of the fluid 
(Q) and the 3-dimensional velocity space (IR 3 ). Moreover, e I are two additional 
real (hidden) variables, with d denoting in particular an ignorable configuration-space 
variable [both for the fluid fields and the kinetic distribution function /(x,£)] defined in 
a bounded interval 1$ = [$oj$i] C K. The streaming operator L in this case reads 
L = i + v-i + w^ + £-{F(x,t)} + £{H(x,t)} 7 where F^x,*) = {F(x, t), ff(x, t)} 
can be interpreted as a mean field force acting on a particle with state x = (r,^, v,w). In the 
sequel we intend to prove that, at least in a suitable finite time-interval J, the fluid fields 
V,p, T can be identified with the velocity moments J rfviG(x,t)/(x,t), where respectively 

G(x,t) = v/po, (v — V) 2 /3, mw 2 /3p and /(x,i) is a properly defined kinetic distribution 

function. In addition, if the same distribution function /(x,t) is strictly positive in the whole 

set T x / and the statistical entropy functional S(f) — — J dx/(x,t) ln/(x,£) exists for all 

r 
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t G I, we intend to show that the thermodynamic entropy can always be identified with 
S(f),i.e., that 

St = S(f). (12) 

To reach the proof, let us first show that, by suitable definition of the "force" fields F(x, t) 
and Q(~x,t), a particular solution of the the IKE (Tj0) is delivered by the (extended-space) 
Maxwellian distribution: 

P I 1 

/ M (x, t) = -5-5 exp < — — > . (13) 

^ V th,p V th,T { V 2 th p vl h T J 

Here u = v — V(r, t) is the relative velocity, while v t h, P = ^/2pi( r , t)/p and v t h,r = 
A/2T(r, t)/m denote respectively the pressure and temperature thermal velocities. Fur- 
thermore pi(r, t) = pq + p(r,t) is the kinetic pressure. In these definitions, Po{t) (to be 
denotes as pseudo-pressure) is an arbitrary strictly positive and suitably smooth function 
defined in /, while the mass m > is an arbitrary real constant. The following theorem can 
immediately be proven: 

Theorem 1 - Local Maxwellian solution for INSFE-IKT Let us identify respec- 
tively the vector and scalar fields F(x,t) and H(x,t) with 

F(x,t;/ M ) = F + F 1 (14) 



H(*,t;fM) = ^K+^u~lnT, (15) 
where Fq, F-^ and K read respectively 

F (x,t; f M ) = -f + ~u • VV+iw • u+z/V 2 V, (16) 
p Q z z 

FiCM; fa) = ^~A + ^Vln Pl - H , (17) 



2pi 2 ^ K 2]' 
_o.„ i/ /<9V* , dV k \ 2 I 



K = X V 2 T+ — (-^ + -^\ +— , (18) 
2c p \ox k dxi J pc p 



while, denoting = + V-V, A reads: 



* D 9 



^-V + V • VV - — f - z/V 2 V 
at p 



(19) 



It follows that: 
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1 ) the local Maxwellian distribution < f73j) is a solution of the IKE if an only if the 
fluid fields {p = p Q > 0, V,p, T} satisfy the fluid equations (GP-j^l); 

2) the velocity-moment equations obtained by taking the weighted velocity integrals of 
Eq. with the weights G(x,t) = 1, v/p a , (v — V) 2 /3, w 2 /3 deliver identically the same 
fluid equations (E]]-([7]). 

PROOF - First we notice that if the fluid equations ([31)-([7| are satisfied identically in Q x J, 
the proof that ffl3|) is a particular solution of the IKE [Eq.(j4])] follows by direct differentiation. 
The converse implication, i.e., the proof that if (Tl3l) is a solution of Eq.(jl]) then the fluid 
equations ©-(El) are satisfied identically in Q x /, follows by evaluating the velocity moments 
of Eq. (1261) for the weights G = 1, v, (v — V) 2 /3, w 2 /3. In analogy to the case of isothermal 
fluids the present theorem can be generalized to suitably smooth non-Maxwellian initial 



distribution function 



3j. In such the fields F , F : and H read respectively 



F (x,t; /) = - [V-S - Vpi + f] + u • VV+z/V 2 V, (20) 

Po 

Fl (x,t; /) = lu \-A+^V ■ Q - \ [V-Hl Q) + ^V-H - 3 -\ , (21) 
while the scalar field if(x, t; /) now reads 

tf(x,t;/) = !^i + fu-!-lnT. (22) 

Here i^i denotes the scalar function K\ = K— ■ V InT + V ■ (V^T) — V • VT] , where ii' 
is defined by Eq. (1181) and furthermore the moments Q, n, $t and Vt are defined respectively 
as Q = / d 3 vdwu ! jf, U = J d 3 vdwuuf, $ T = / d^vdw^uf and V T = J d^vdw^wf. 



III. H -THEOREM AND THE 2ND PRINCIPLE OF THERMODYNAMICS 

In this section we want to prove that, provided the pseudo-pressure po(t) is suitably 
defined (i.e., is a uniquely-prescribed function of time), an H Theorem can be established 
for the statistical entropy S (/) = — j T dx.f In /, which warrants the strict positivity of the 
kinetic distribution function in the whole set T x /. As a further consequence, the position 
(j!2p holds too. For this purpose we distinguish between isothermal and non-isothermal 
fluids, i.e., fluids in which the fluid temperature is respectively a constant in the whole 
set Q x J, or not. In particular, we intend to prove that, in the first case (isothermal 
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fluid), a constant H-theorem holds for the statistical entropy under suitable conditions, i.e., 
by suitably prescribing the pseudo-pressure. Instead, to reach an H-theorem for a non- 
isothermal fluid, it is necessary to include also a suitable prescription on the r.h.s. of the 
Fiourier equation [and in particular on the scalar field J; see Eq.(j7J)], which defines the 
quantity of heat generated by external sources. In both cases, the result can be proven to 
hold at least in a finite time interval / and for an arbitrary strictly positive (and suitably 
summable) distribution function /. 

Theorem 2 - H-theorem Let us assume that: 1) O is a bounded subset of M 3 ; 2) 
the kinetic distribution function f coincides identically in Q x I with the local Maxwellian 
distribution f m defined by Eg. fTS)) . Furthermore, let us distinguish respectively the cases in 
which the fluid is isothermal in Q x I or not. In the first case we demand that the following 
assumption is fulfilled: 3 A) the pseudo-pressure po(t) is po(t) > determined in such a 
way to satisfy identically Vt E I the constraint 



dr — 

n Pi 



d 1 

^-pi + V-Q Vp-Q 

at pi 



0. (23) 



Instead, for non-isothermal fluids, we require that the following two assumptions are satisfied 
(SB and 4)- 3B) the pseudo-pressure Po(t) is Po(t) > and satisfies identically Vt G / the 
constraint 



3 (a 1 



JU + V-Q--Vp-Q 

at pi 



2 Jn Pi 
l - Po J rfri [$ T • V In T + V ■ (V T T)] = 



(24) 



4 ) the quantity of heat generated by external sources J, either vanishes identically in Q x I 
(isolated fluid) or Vi G / is externally heated in the sense of the inequality . 

Then it follows respectively: A ) for isothermal fluids: the statistical entropy S (f) is 
constant, i.e., there holds identically the constant H-theorem: 

J^(/) = 0; (25) 

B) for non-isothermal fluids: the statistical entropy S (f) is a monotonic function of time, 
i.e., it holds, instead, the H-theorem\/t G / 

|s(/)>0. (26) 
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PROOF - For definiteness let us first consider an isothermal fluid, i.e., requiring T = const. 
in Q x I. In this case the proof of Eq. fl25l) is immediate. In fact for an arbitrary, suitably 
smooth and non vanishing, non-Maxwellian distribution the entropy production rate j^S (/) 
results 



3 I A 1 

~Po / ar- 



^Pi + V • Q--Vj?-Q 
In Pi lot p 

which thanks to the constraint equation Eq. (l23|) for po(t) vanishes identically. In particular, 
in the case in which the kinetic distribution function / coincides with the local Maxwellian 
distribution (J4]) the same equation delivers 



d 3 
Wt S{f) = 2 P < 



f d 

/ dr-\n[p (t)+p(r,t)] = 0. (28) 



Let us now consider the case of a non-isothermal fluid obeying the fluid equations ([5])- (JTJ) . 
In this case, imposing onpo(^) the constraint (!24l) the entropy production rate reads 

yW-\f.j*K. (29) 

Hence, for a fluid which is either isolated or subject to external heating, in the sense of the 
inequality (Q, the H-theorem (j26|) manifestly holds. 

This result is consistent with the second principle if thermodynamics [i.e., the inequality 
(JSJ)]. As a consequence, this enables us to specify also the thermodynamic entropy in terms 
of a suitable phase-space moment of the kinetic distribution function /. 



IV. CONCLUDING REMARKS 



A basic implication of the IKT here developed for INSFE is that it has been constructed 
in such a way to satisfy the following requirements: 

1. completeness: all fluid fields are expressed as moments of the kinetic distribution func- 
tion and all hydrodynamic equations can be identified with suitable moment equations 
of IKE; 

2. closure condition of moment equations: there must exist a subset of moments of IKE 
which form a complete system of equations, to be identified with the prescribed set of 
hydrodynamic equations; 
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3. smoothness for the fluid fields: the fluid fields are assumed suitably smooth so that the 
solution of the kinetic distribution function exists everywhere in a suitable phase-space; 

4. arbitrary initial and boundary conditions for the fluid fields: the initial conditions for 
the fluid equations are set arbitrarily while Dirichlet boundary conditions are consid- 
ered on the boundary; 

5. self- consistency: the kinetic theory holds for arbitrary (and suitably smooth) initial 
conditions for the kinetic distribution function. 

6. non- asymptotic IKE: i.e., the correct hydrodynamic equations must be recovered by 
the inverse kinetic theory independently of any dimensionless parameter. 

The present approach has the following main features: 

1. a suitable classical dynamical system has been constructed which uniquely determines 
the evolution of the fluid fields; 

2. the IKT is based on an extended phase-space formulation which relies of the micro- 
scopic statistical description of the dynamical system; 

3. the theory satisfies the second principle of thermodynamics. 

An interesting result of the theory, relevant for the mathematical investigation of the fluid 
equations, concerns the discovery of the underlying dynamical system, i.e., the phase-space 
classical dynamical system ffTUl) . We have found that this can be identified with a - gen- 
erally non-conservative - dynamical system advancing in time the microscopic distribution 
function and generated by the kinetic equation itself. The evolution of the fluid fields is thus 
determined uniquely by this dynamical system, a result that in principle may be achieved 
without solving explicitly the fluid equations themselves. 

Acknowledgments 

Work developed in cooperation with the CMFD Team, Consortium for Magneto-fluid- 
dynamics (Trieste University, Trieste, Italy). Research developed in the framework of the 
MIUR (Italian Ministry of University and Research) PRIN Programme: Modelli della teoria 



10 



cinetica matematica nello studio dei sistemi complessi nelle scienze applicate. The support 
of COST Action P17 (EPM, Electromagnetic Processing of Materials) and GNFM (Na- 
tional Group of Mathematical Physics) of INDAM (Italian National Institute for Advanced 
Mathematics) is acknowledged. 

Notice 

§ contributed paper at RGD26 (Kyoto, Japan, July 2008). 



11 



[1] M. Ellero and M. Tessarotto, Bull. Am Phys. Soc. 45 (9), 40 (2000). 

[2] M. Tessarotto and M. Ellero,Proc. 24th RGD, Bari, Italy (July 2004), Ed. M. Capitelli, AIP 

Conf. Proc. 762, 108 (2005). 
[3] M. Ellero and M. Tessarotto, Physica A 355, 233 (2005). 

[4] M. Tessarotto and M. Ellero, Proc. 25th RGD (International Symposium on Rarefied gas 
Dynamics, St. Petersburg, Russia, July 21-28, 2006), Ed. M.S. Ivanov and A.K. Rebrov 
(Novosibirsk Publ. House of the Siberian Branch of the Russian Academy df Sciences), p. 1001; 



arXiv:physics/0611113 (2007). 



[5] M. Tessarotto and M. Ellero, Physica A 373, 142 (2007); arXiv: [physics/0602140 



[6] M. Tessarotto, M. Ellero and P. Nicolini, Phys. Rev. A 75, 012105; arXiv:quantum-ph/060691 
(2007). 

[7] M. Tessarotto, M. Ellero, M. Mond, N. Asian and P. Nicolini, An exact pressure evolution 



equation for the incompressible Navier-Stokes equations, |arXiv:physics/0612072| contributed 
paper at RGD26 (Kyoto, Japan, July 21-25, 2008). 
[8] S. Succi, The Lattice- Boltzmann Equation for Fluid Dynamics and Beyond (Numerical Math- 
ematics and Scientific Computation), Oxford Science Publications (2001). 
[9] C. Cremaschini, A. Beklemishev, J.C. Miller and M. Tessarotto, Axisymmetric gravitational 
MHD equilibria in the presence of plasma rotation, contributed paper at RGD26 (Kyoto, 
Japan, July 21-25, 2008). 

[10] M. Tessarotto, E. Fonda and M. Tessarotto, The computational complexity of traditional 
Lattice- Boltzmann methods for incompressible fluids, contributed paper at RGD26 (Kyoto, 
Japan, July 21-25, 2008). 

[11] E. Fonda, M. Tessarotto and M. Ellero, Lattice- Boltzmann inverse kinetic approach for clas- 
sical incompressible fluids, contributed paper at RGD26 (Kyoto, Japan, July 21-25, 2008). 

[12] Marco Tessarotto, Claudio Cremaschini, Piero Nicolini and Massimo Tessarotto, Lagrangian 
dynamics of incompressible thermofluids, contributed paper at RGD26 (Kyoto, Japan, July 
21-25, 2008). 



12 



